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We present a study of electric transport at high temperature in a model of strongly interacting 
spinless fermions without disorder. We use exact diagonalization to study the statistics of the energy 
eigenvalues, eigenstates, and the matrix elements of the current. These suggest that our nonrandom 
Hamiltonian behaves like a member of a certain ensemble of Gaussian random matrices. We calculate 
the conductivity o(u>) and examine its behavior, both in finite size samples and as extrapolated to 
the thermodynamic limit. We find that o-{uj) has a prominent non-divergent singularity at u = 
reflecting a power-law long-time tail in the current autocorrelation function that arises from 
nonlinear couplings between the long-wavelength diffusive modes of the energy and particle number. 
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I. INTRODUCTION 

The subject of transport in strongly-correlated electron 
systems is of considerable current interest, motivated, for 
example, by recent experiments on cuprate and cobal- 
tate materials measuring electrical and thermal transport 
properties P, . Traditionally, much of the work has fo- 
cused on the behavior of quantum systems at low tem- 
peratures Q, where transport properties can be used to 
characterize a material's ground state. Under these con- 
ditions, the excitations responsible for transport are of- 
ten dilute and weakly interacting ( "quasiparticles" ) , thus 
allowing for a straightforward description of transport, 
e.g. using the Boltzmann (kinetic) equation formalism 
[J] . Interestingly, in the opposite extreme, when temper- 
ature is high and no well denned elementary excitations 
exist, transport properties remain nontrivial and diffi- 
cult to calculate and understand [j| . Additionally, what 
makes studying this regime worthwhile is the fact that 
sometimes interesting phenomena take place over a broad 
temperature range, extending up to temperatures so high 
that direct interpretations in terms of low temperature 
universal properties (either conventional or exotic) can 
be questioned, thus leaving the high temperature regime 
addressed in this work as a promising starting point of 
analysis. Although we attempt no detailed comparisons 
with experiments here, it may be of interest to note that 
the optical conductivity we find is very broad, extend- 
ing over frequencies of order the bare bandwidth, with 
DC resistivity growing linearly with temperature. These 
features are remarkably reminiscent of what is often ob- 
served in the normal state of correlated materials as di- 
verse as organic salts, CDW systems, Ceo and even high 
temperature superconductors. However, further work is 
necessary to elucidate possible implications of our high 
temperature approach to any specific material. 
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In this paper, we study the conductivity, cr, as a func- 
tion of frequency to and temperature T. For a finite quan- 
tum system with a discrete spectrum, the real part of the 
conductivity at nonzero frequency is given by the Kubo 
formula: 
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a{uj,T) = TT Y j e- flE -\J nm \ 2 5{E n -E m -nw) , 
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(i) 

where n and m are eigenstates of the Hamiltonian, E n 
and E m are the corresponding energy eigenvalues, J nm is 
the matrix element of the total current operator between 
these two states, (3 — l/fc^T and Z is the partition func- 
tion. Note that we do not include any "external" dissi- 
pative bath; what we are studying is the absorption of 
energy from an infinitesimal applied AC electric field by 
transitions between eigenstates of the Hamiltonian. 

One thing we investigate here is how the above singu- 
lar expression for the conductivity converges to a contin- 
uous function in the thermodynamic limit. It is a sum 
of delta- functions, one for each pair of energy eigenstates 
that are connected by the current operator. The num- 
ber of eigenstates scales exponentially in the number of 
degrees of freedom in the system ~ e sL , where L is sys- 
tem's linear extent and s is the entropy density, with 
average level spacing A ~ e~ sL . For a generic non- 
integrable Hamiltonian there are only a few good quan- 
tum numbers. These include the total number of par- 
ticles, the Bloch (lattice) momentum for systems with 
discrete translational invariance, and in some models of 
interest the total spin. All pairs of states with the same 
quantum numbers generically contribute a delta-function 
to the conductivity, and so the total number of terms in 
the sum scales as ~ e 2sL . Thus the number of delta- 
functions in a(u>) grows very rapidly both with increasing 
system size and with increasing temperature (and thus 
s), and away from zero frequency a good approximation 
to the thermodynamic limit is rapidly approached with 
only very fine binning or other smoothing of the delta- 
functions. However, there are significant and interesting 
finite size effects present near zero frequency, and eluci- 
dating them is one of the main topics of this paper. 
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For models where a finite-sized system has a bounded 
spectrum and a finite-dimensional Hilbcrt space, such as 
the one-band model of interacting lattice fermions that 
we will study here as our example Hamiltonian, the above 
expression simplifies in the high temperature limit to 

k B Ta(cj) \ J nm\ 2 S(E n - E m - tiu) . (2) 

We will mostly focus on this limiting behavior of the 
conductivity, since it is already nontrivial. However, all 
of the features we will discuss appear to apply at finite 
temperatures as well. One basic noteworthy feature of 
this high T limit is that the resistivity is proportional to 
the temperature. 

Such a simple statistical characterization of the con- 
ductivity breaks down at low temperatures and frequen- 
cies. At low temperatures both thermodynamics and dy- 
namics are dominated by comparatively few transitions 
from low energy states residing in the tail of the Hamilto- 
nian's full spectrum, and their properties may depend on 
specific details of the system and its ground state. The 
low frequency part of the high temperature conductiv- 
ity, on the other hand, shows various finite size effects 
that can be understood based on rather general consid- 
erations: 

i) Level repulsion reduces the number of pairs of states 
with a very small energy difference, producing an under- 
estimate of the macroscopic low-frequency conductivity 
(in fact, a vanishes as w — > 0); 

ii) If the above formula (f) for the conductivity of a 
finite system is applied precisely at u = 0, we find a delta 
function due to the current's diagonal matrix elements 
J nn . For integrable Hamiltonians this delta- function has 
a nonzero weight in the thermodynamic limit, hence the 
macroscopic DC conductivity appears to be infinite at 
all temperatures. In weakly non-integrable systems this 
delta function is weakly broadened, it is the celebrated 
Drude peak. In our case, and in non-integrable systems 
more generally, the weight of this precisely elastic delta- 
function in the conductivity scales exponentially to zero 
with increasing system size, as discussed in Section IV 
below. 

iii) Finally, in the infinite system there is a non- 
divergent zero- frequency singularity of the form a(uj — 
0) — cr(u>) ~ oj d / 2 due to nonlinear couplings between 
the long- wavelength diffusive energy and particle density 
modes. This singularity is rounded by finite-size effects. 

We discuss and disentangle all of these effects (as 
well as demonstrate the general validity of considerations 
away from u) = 0) using a simple one dimensional model 
as a case study. While some of the results are special to 
the particular model, e.g. the dependence of the low fre- 
quency non-analyticity on dimension and range of inter- 
actions, the overall picture that emerges is clear and ro- 
bust against variations in parameters, and most features 
are not special to one dimension. It also appears to ap- 
ply at all temperatures in the disordered phase provided 



/3a; « 1. In Section II we present the Hamiltonian un- 
der study as well as the numerical procedures employed 
to solve the problem exactly for finite sizes. In Section 
III we use the exact numerical results to study statistical 
properties of the spectrum and the eigenstates, including 
the matrix elements of the current operator. The con- 
ductivity is computed and discussed in Section IV, and 
its diffusion-induced non-analyticity is analyzed in Sec- 
tion V. Finally, we conclude with a brief overview and 
list of open questions. 

II. t - t' - V MODEL AND EXACT 
DIAGONALIZATION 

Although much of what we find is rather generic to 
non-integrable quantum many-body solid-state systems, 
the specific model we consider here is a one dimensional 
chain of spinless fermions with a single tight-binding or- 
bital at each site. The Hamiltonian is 

3 3 3 

(3) 

where dj, Cj and rij = CjCj are fermion creation, anni- 
hilation and density operators, respectively. This is one 
of the simplest non-integrable (quantum chaotic) many- 
particle models. For spinless lattice fermions the size of 
the Hilbert space is small, with only 2 states per site, 
and thus 2 L states for a chain of length L sites. We have 
a short-range (nearest-neighbor) interaction V (spinless 
lattice fermions cannot have an on-site interaction). The 
second-neighbor hopping t' makes the Hamiltonian non- 
integrable and also breaks particle-hole symmetry. This 
system's good quantum numbers (conserved quantities 
in addition to the energy) are total particle number and 
total crystal momentum. As written above, with the pe- 
riodic boundary conditions that we use, it also has parity 
symmetry. We have chosen to break parity with an irra- 
tional phase twist at the boundary (equivalent to thread- 
ing a fraction of a magnetic flux quantum through the 
loop). This minimizes the number of distinct symmetry 
sectors to those labeled by the two good quantum num- 
bers. We have also looked at the model where the second- 
neighbor hopping is replaced by a second-neighbor inter- 
action; this latter model is equivalent to that studied by 
Rabson, et al. Q}. We chose to focus instead on our 
model because we did not want particle-hole symmetry. 
However, the results we present are generally indepen- 
dent of the presence or absence of any of these sym- 
metries other than conserved particle number. Similar 
models have been used in earlier studies of charge, heat 
and/or spin transport and the crossover from integrabil- 
ity to non-integrability |a, |6( . 

For this model the operator for the total particle cur- 
rent is 

3 = | Yl c t c 3+i + ^J-t' E S +c J+2 + h.c. (4) 

3 3 
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We have also looked at heat transport and thermopower, 
with similar preliminary results; this work will be re- 
ported later. We work for convenience at half filling, and 
average over the different total crystal momentum sec- 
tors. It appears, as we expect, that the behavior is sta- 
tistically the same in each such momentum sector, and 
that there is nothing special about half filling for the 
regimes and properties we examine. We specifically take 
t = t' = 1 and V = 2, so the interaction is roughly equal 
to the single-particle bandwidth, and thus strong. All 
energies are measured in units of t = t' = 1, frequencies 
in units of t/h — 1, and lengths in units of the lattice 
spacing. 

We use a Householder diagonalization routine to ex- 
actly diagonalize the Hamiltonian within each momen- 
tum sector (typically 16 or 18 sites, with 8 or 9 parti- 
cles, respectively). For convenience, we choose to work in 
the basis of single-particle momentum eigenstates. These 
states are also eigenstates of the total current operator, 
so this is one of the two "natural" bases for this prob- 
lem. The other natural basis is less simple; it is that of 
the eigenstates of the other operator that enters in deter- 
mining the conductivity, namely the Hamiltonian itself. 



III. STATISTICS OF THE SPECTRUM, 
EIGENSTATES AND CURRENT MATRIX 
ELEMENTS 

One of the common measures used to characterize the 
spectrum of a non-integrable quantum Hamiltonian is 
the (Wigner-Dyson) distribution of level spacings. Our 
Hamiltonian is a real matrix (this is true in our basis of 
single-particle momentum eigenstates even with a mag- 
netic flux through the loop), so the probability distri- 
bution of the level spacing should be that of the GOE 
ensemble of random matrices, 

P(s ) = 2 /4 ; (5) 

where s here is the level spacing in units of the mean level 
spacing at that energy. To compare to this prediction we 
first have to measure the density of states, which sets the 
mean level spacing. The average density of states N(E) 
in each sector of total crystal momentum is shown in Fig. 
1 for the case of chain length L = 18. Thus we obtain an 
energy-dependent mean level spacing A(E) = l/N(E). 
The spacing of each pair of adjacent levels within one 
sector is divided by A(E) to yield the scaled level spacing 
s. The resulting probability distribution P(s) for L = 18 
is shown in Fig. 1, where we can see it agrees very well 
with the GOE form ©. 

The Kubo formula shows that a relevant statisti- 
cal quantity that enters the conductivity at a given fre- 
quency uj is the mean-square value of the current matrix 
elements J mn between energy eigenstates separated by 
an energy hui = \E m — E n \. [Note that our eigenstates 
all have real amplitudes, so these matrix elements are all 
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FIG. 1: (Top) The average density of states in each symmetry 
sector. (Bottom) The probability distribution of the scaled 
level spacing for the L = 18 system (data points) and the 
GOE form for this distribution (line), showing the excellent 
agreement. 



real here.] Thus we were interested in characterizing the 
probability distribution of these current matrix elements. 
What we find is that the distribution is consistent with 
a Gaussian, with mean zero and a variance that depends 
on both the initial and final energies. The variance has a 
strong dependence on ui, giving <j(lo) its dependence on ui 
(see below). The variance also has a weaker dependence 
on the average energy E = (E m + E n )/2. To examine 
the probability distribution of J mn , we bin the matrix 
elements by both u> and E, and look at all the J mri 's 
that fall in one such bin of energies, accumulated over all 
total momentum sectors. We measure the first through 
fourth moments, Mi through M4, of these distributions, 
observing that the odd moments are consistent with zero, 
and looking particularly at the ratio R = M4/M2 that 
is equal to 3 for the Gaussian distribution. If the bins 
are chosen too wide, the variance changes over the bin 
and the resulting distribution has a ratio R larger than 
3. But as the bins are made narrower, R converges to 
3, indicating that the current matrix elements do indeed 
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have Gaussian probability distributions. Since the vari- 
ance depends much more strongly on to than on E, the 
bins are taken to be much narrower along the u> direction. 

Thus we conclude that the total current operator, as 
written in the basis of the eigenstates of the Hamilto- 
nian, is effectively a Gaussian random matrix that is, 
loosely speaking, "banded" , meaning the variance of the 
matrix elements J mn depends smoothly on E m and E n . 
The "bands" of constant variance in this case do not 
run exactly parallel to the diagonal of the matrix. Of 
course, neither our Hamiltonian nor the current operator 
contain any random variables; the apparent randomness 
comes from the statistical properties of the highly ex- 
cited eigenstates of our nonintegrable (and thus "quan- 
tum chaotic" ) Hamiltonian. The total current operator, 
on the other hand, is integrable, with its eigenstates be- 
ing the simultaneous eigenstates of all the single-particle 
momenta. When we write the Hamiltonian as a matrix 
in the basis of the eigenstates of the current operator, 
which is precisely what we do to diagonalize it numeri- 
cally, it is a highly regular sparse matrix. [We have also 
looked at the energy current operator, which enters in 
the thermal conductivity and thermopower; the energy 
current operator is nonintegrable, and the Hamiltonian 
does look like a Gaussian random matrix when written 
in terms of the eigenstates of the total energy current.] 

The final statistical characterization we have done is 
of the amplitudes (n\a) when the eigenstates \n) of H 
are written in the basis of the eigenstates |a) of J (and 
vice versa). These amplitudes, like the current matrix 
elements J mn discussed above, appear to be mean-zero 
Gaussian random variables with a variance depending 
on two energies, namely E n and the kinetic energy K a 
of the current eigenstate. Each current eigenstate is also 
an eigenstate of the kinetic energy operator, which is the 
sum of the hopping terms in the Hamiltonian. As ex- 
pected, a state \n) of high (low) total energy E n consists 
primarily of states \a) of high (low) kinetic energy K a , 
but it appears to be quite uniformly and randomly ex- 
tended over all the current eigenstates at each particular 
kinetic energy. Thus we find that both the eigenstates of 
H , and their current matrix elements appear to behave as 
Gaussian random variables with energy-dependent vari- 
ances. 
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FIG. 2: (Top) Drude weight as a function of L. The Drude 
weight is expected to vanish exponentially with L. (Bottom) 
Conductivity as a function of frequency for L = 18, L = 16 
and L = 14. The agreement among different system sizes 
is good at high frequencies but significant finite-size effects 
can be seen at low frequency. (Inset) Closer view of the low 
frequency regime. The dip in the conductivity around u = 
is due to level repulsion. This dip rapidly narrows with 
increasing system size. The bins used here have a width equal 
to the mean level spacing at the maximum of the density of 
states for that L, except in the case of the main panel for 
L = 14, where the bins widths are instead chosen proportional 
to ui in order to reduce the scatter. 



IV. THE CONDUCTIVITY 

The conductivity is given by the Kubo formula Q. 
As remarked earlier, it is a series of delta functions for 
a finite-sized system, which can be binned to produce a 
smooth curve for <j(lo). Since the limit of zero frequency 
is of particular interest, we will focus much of our at- 
tention on low frequencies. Although the expression (1) 
does not apply at strictly zero frequency, it does con- 
tain a delta- function there that all levels contribute to. 
This zero-frequency feature has often been called (am- 
biguously) the "Drude peak" , and appears to have finite 



weight in the thermodynamic limit for integrable mod- 
els. However, for nonintegrable quantum chaotic models 
such as we study here, the weight in this strictly elas- 
tic Drude peak vanishes strongly in the thermodynamic 
limit, see Fig. 2, top panel. This can be roughly un- 
derstood by noting that the full conductivity (1), which 
has a finite large L limit, is a sum of ~ e 2sL delta func- 
tions, where s here is the entropy density. Of these, only 
e sL are at zero frequency (since there are no degeneracies 
within any symmetry sector). We find that the weight of 
the zero frequency terms in the sum are of the same or- 
der as those at small frequency, so the relative strength 
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of the so-called Drude peak at zero frequency drops as 
e~ sL with increasing L. A more careful analysis gives 
power-law in L prefactors to these exponential depen- 
dences. But the presentation on a semilog plot in Fig. 2 
shows that our results are consistent with a Drude weight 
that vanishes exponentially with increasing L. 

The conductivity is shown in Fig. 2 for our longer 
chains, of length 14, 16 and 18. To obtain fairly smooth 
functions, the delta functions are binned (see caption), 
and the elastic uj — "Drude" peak is left out of these 
plots. At high frequencies, the finite size effect is clearly 
very small, so these short chains give a good approxima- 
tion to the thermodynamic limit. An interesting question 
that we do not yet have an answer for is: What is func- 
tional form of the conductivity at high frequencies (say, 
u) > t), and what is the physics determining it? We find 
that <j{uj) appears to decrease faster than any power of 
lu in this regime. 

While there appears to be good convergence to the 
thermodynamic limit of cr(w) at high frequencies for these 
chain lengths, clear finite-size effects are seen at low fre- 
quencies. We find two sources of these finite-size effects: 
level repulsion, and diffusive "long-time tails" . The for- 
mer is expected and has been discussed before in this and 
many other contexts, but the latter appears to be a new 
observation for this type of quantum many-body system. 
Level repulsion causes a reduction of <t(w) at frequencies 
of order or less than the mean level spacing, because there 
are fewer pairs of levels with those energy differences that 
can produce the dissipation. The level spacing vanishes 
exponentially in the system size, so this feature, which 
appears as a dip in <j(uj) around zero frequency, narrows 
very rapidly with increasing L, as is apparent in the inset 
in Fig. 2. 

Naively, one might expect (as we did when we started 
this project) that in this high temperature disordered 
regime with no long-range spatial correlations the con- 
ductivity of the infinite system would be a smooth an- 
alytic function of frequency in the vicinity of zero fre- 
quency. However, examination of our results in Fig. 2 
show that this does not appear to be the case. The con- 
ductivity appears to have a sharp maximum at u> = 
that gets steadily sharper as L is increased. This sharp 
peak is apparent in published numerical data from other 
one-dimensional quantum models but the reason for 
it appears not to have been discussed. As we show in the 
next section, long-wavelength diffusive modes of the en- 
ergy and particle density generically interact nonlinearly 
to produce a nondivergent zero-frequency singularity of 
the form <j{lo) = a — &ypJ] + ... in these one-dimensional 
systems. In the top panel of Fig. 3 we show that our 
data appear to be consistent with such a form for the 
low- frequency conductivity in the large L limit. This 
behavior arises from a "long-time tail" in the autocorre- 
lation function of the total current that decays with time 
difference as (J(O)J(t)) ~ i~ 3 / 2 at long times. We show 
in the bottom of Fig. 3 that our data for this autocorre- 
lation for L = 18 do show such a power-law decay over a 
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FIG. 3: (Top) Conductivity plotted versus ^/uj for L = 16 
and L = 18. The straight line is a linear fit on this plot. The 
conductivity appears linear in ^/uj at low frequencies before 
it rounds off at very low frequency due to finite-size effects. 
The finite-size effects are, as expected, more pronounced for 
L = 16 than L — 18. (Bottom) Autocorrelation function 
of the total current vs. time on a log-log plot. There is a 
long-time tail with (J(O)J(t)) ~ t~' i/2 . 



significant time range. 



V. LONG-TIME TAILS DUE TO DIFFUSIVE 
MODES 

Here we first give a more general argument inspired 
by our one-dimensional system. Consider an infinite 
many-body system in a disordered phase. If it has any 
conserved densities, such as energy, momentum, angu- 
lar momentum, particle number etc., we may ask about 
the transport properties for these quantities. A long- 
wavelength disturbance in the conserved quantities re- 
laxes at a (possibly complex) rate T(q) where q is the 
wavevector of the disturbance, and this rate is in general 
different for different linear combinations of the various 
conserved quantities. In the absence of long-range inter- 
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actions or propagating modes, the relaxation is diffusive: 
r ~ Dq 2 for small q. These long-wavelength diffusive 
modes in general produce singular behavior of the trans- 
port properties in the zero- frequency (long-time) limit, as 
discussed, for example, by Kirkpatrick, et al. Q. This 
has been discussed in the context of fluids where mo- 
mentum conservation makes these effects quite strong, 
and in systems where quenched randomness appears to 
play an essential role 0. 

For the Hamiltonian that we are studying, there are 
two conservation laws that give diffusive modes: energy 
and particle number. The lattice (umklapp) breaks the 
momentum conservation down to just the conservation 
of crystal momentum, which does not appear to produce 
any slow modes in the disordered phase. 

To be a little more general, let us first consider a 
translationally- and rotationally-invariant system in a 
disordered phase with some number of conserved den- 
sities n a (r,t) in possibly more than one dimension. The 
conservation laws then dictate 



where j Q (r, t) is the current of conserved quantity a. Let 
us define the local densities the deviations from 

the average densities, so that our system has zero aver- 
age densities (n a ) — 0. We assume the long- wavelength, 
low- frequency dynamics is diffusive, and expand the equi- 
librium coarse-grained dynamics in the densities and the 
wavevectors (spatial gradients). The currents are ex- 
panded as 

ia = -D al3 S7np + Ea^npVn^ + ... + T] a , (7) 

where D is the diffusivity matrix, r/ is noise (due to non- 
linear couplings to short- wavelength modes) with an in- 
tensity given by a fluctuation-dissipation relation (e.g., 
0), and repeated subscripts are summed over. The 
lowest-order (in density and gradient) nonlinear correc- 
tion to the diffusive behavior has also been shown explic- 
itly; for example the term E a py corresponds to the linear 
dependence of the local diffusivity matrix element D ai 
on the local density np. 

What should be the low-frequency behavior of a trans- 
port coefficient like <t(uj)? By the Kubo formula, this is 
given by the time fourier transform of the autocorrelation 
function of the total (zero wavevector) current. Integrat- 
ing the above expression, and thus leaving out any terms 
that are total derivatives, the total current is 

Jqi(*) = / iadr = 
f[ E <*»- E ~"' (ngVn 7 - i^Vng) + ... + r) a ]dr . (8) 

Note that the total current is correlated with the densi- 
ties only through the nonlinear corrections to simple dif- 
fusion; we will consider the generic case when these cor- 
rections are indeed nonzero. Through this nonlinearity, 
the zero- wavevector current is coupled to the slow, small- 
wavevector diffusive density modes; this is what produces 



the long time tails. However, since the leading nonlinear 
coupling that we show explicitly here must be to an an- 
tisymmetric combination of the conserved densities, for 
a system with only one conserved density (and thus no 
antisymmetric combinations of densities), one has to go 
to significantly higher-order nonlinearities (|Vn| 2 Vn) to 
obtain the leading long-time tail in the current autocor- 
relations. 

Writing the total current in terms of the long- 
wavelength density fluctuations in reciprocal space, we 
have 

Ja = Ea(il 2i EalP J <V r i l3 ((i)n 7 (-(i)dci + ... + r] a , (9) 

where H/a(q) and n 7 (q) are the spatial Fourier transforms 
of n ) 3(r)and n 7 (r). The low momentum fluctuations in 
n ( g(q) and n 7 (q) decay diffusively. The eigenmodes of 
that decay are the linear combinations of the densities 
that diagonalize the diffusivity matrix. What we are do- 
ing here is considering the correction due to nonlinearity 
to the linear diffusive "fixed point" theory. These non- 
linearities are irrelevant in a renormalization group ap- 
proach, but give singular corrections to the leading dif- 
fusive behavior. The long-time tail we discuss here is 
one of those corrections. At this diffusive fixed point, in 
the basis of densities that diagonalizes the diffusivity ma- 
trix, the autocorrelations of the long-wavelength density 
fluctuations decay as 

(n / j(q,0)n 7 (q / ,t)> ~ e-^'fy^q + q') (10) 

for small q and large t, where Dp is the diffusivity for 
eigenmode f3. Thus at long times, the autocorrelation of 
the total current behaves as 

(J(O)J(t)) - J d(\q 2 exp{-Dq 2 t) + ... 

~ l/t^ 1 )/ 2 (11) 

in d dimensions, where D = Dp+D 1 is the sum of the two 
smallest diffusivities. Fourier transforming this to get the 
conductivity, we see that for one dimension and small uj, 
cr(u>) ~ a — by/\uj\ + ...; our numerical results (Fig. 3) are 
quite consistent with this. In our finite-sized systems, the 
integrals on q are replaced by sums, which has the effect 
of rounding out this singularity at a frequency of order 
D/L 2 . This rounding is apparent in Fig. 3, and it can 
be seen to be reduced with increasing L. Note that the 
characteristic frequency of this diffusive finite-size effect 
vanishes only as L~ 2 with increasing size, in contrast to 
the scale for the level repulsion, which vanishes exponen- 
tially in L. 

This explanation of the long-time tail in the autocor- 
relation of the total current makes certain assumptions 
and predictions about the dynamics and correlations of 
the long wavelength modes, which we have checked nu- 
merically within our model. First, it assumes that the 
dynamics of the long-wavelength fluctuations in the par- 
ticle density n and the energy density e are diffusive. We 
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have checked this by calculating their dynamic correla- 
tions: (|n(<7,cj)| 2 ), (|e(g,w)| 2 ), and (n(q, u)e(-q, -u>)) for 
the chain of length L — 18. At small q and ui these fit 
quite well to the expected diffusive behavior. The two 
diffusivities differ by approximately a factor of two. The 
faster eigenmode of the diffusion has a larger component 
of particle density than energy density, and vice versa for 
the slower mode. 

Another assumption is that the long-wavelength modes 
interact nonlinearly to contribute to the total (zero mo- 
mentum) current, as given by Eq. 10 above. For our 
finite-size Id system with discrete momenta and just the 
two densities, n and e, this may be rewritten for the par- 
ticle current as 

J = C^i{svaq){n{q)e{-q)-n{-q)e{q)) + ... + r ] , (12) 

q>Q 

where C is the nonlinear coupling. Thus we have mea- 
sured C to check that it is indeed nonzero. To do this, 
we define the nonlinear combination of the densities at 
momentum q: 

N{q) = i(sinq){n(q)e(-q) - n(-q)e(q)) , (13) 

and measure its correlations with itself (N(q)N(q')) and 
with the current. We find that the correlations between 
iV's at different g's are quite small compared to the vari- 
ances (N 2 (q)). To get estimates of the coupling C, we 
thus note that 

(N(q)J) - C{N*(q)) ; (14) 

this gives us an estimate of C for each q, and we find 
that these estimates agree well with one another at small 
q and are indeed nonzero, as expected. 

The specific form of the long-time tails that we derive 
above are for the case of relaxational dynamics and short- 
range interactions, where the conserved densities relax 
diffusively. Of course, charge carriers in real materials 
interact via the long-range Coulomb interaction, and the 
long-wavelength mode corresponding to the charge den- 
sity, the plasmon, does not show diffusive relaxation; it is 
likely overdamped at high T, relaxing at a nonzero rate. 
Thus conserved energy and charge are not alone suffi- 
cient to produce the long-time tails discussed above once 



Coulomb interactions (or other sufficiently long-range in- 
teractions) are included. What is needed is at least two 
conserved densities that relax diffusively. One possibility 
is ionic conductors, where there are multiple conserved 
charge carriers that can combine to produce another neu- 
tral, diffusively-relaxing hydrodynamic mode in addition 
to the energy. 



VI. CONCLUSIONS 

We have studied the frequency-dependent conductiv- 
ity of a non-integrable fermionic chain with short range 
interactions. Our studies show that such a system can be 
regarded as quantum chaotic in the sense of random ma- 
trix theory, as can the current matrix elements that go 
into the Kubo formula for the conductivity. We have 
confirmed that charge transport in this system is not 
ballistic but diffusive with non-linearities. These non- 
linearities give rise to a non-analytic behavior for the op- 
tical conductivity at low frequencies which we understand 
through hydrodynamic arguments. We observe that the 
system shows finite-size effects at low frequency arising 
from Wigner-Dyson level repulsion and from the low mo- 
mentum fluctuations of the conserved quantities in the 
Hamiltonian. 

The methods used to obtain these results are readily 
generalizable to further explorations of the finite tem- 
perature behavior. The list of open questions includes 
other dynamic response functions and transport coeffi- 
cients (e.g. Peltier and heat conductivities), an extension 
to two dimensions, a systematic investigation of the high 
temperature series for transport, whose first term for the 
conductivity has been computed here, and an extension 
to long range interactions. We shall report on some of 
these in the near future. 
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